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In this paper we study the origin of inertia in a curved spacetime, particularly the spa- 
tially flat, open and closed Friedmann universes. This is done using Sciama's law of 
inertial induction, which is based on Mach's principle, and expresses the analogy be- 
tween the retarded far fields of electrodynamics and those of gravitation. After obtaining 
covariant expressions for electromagnetic fields due to an accelerating point charge in 
Friedmann models, we adopt Sciama's law to obtain the inertial force on an accelerating 
mass m by integrating over the contributions from all the matter in the universe. The 
resulting inertial force has the form F = —kma, where k < 1 depends on the choice of 
the cosmological parameters such as Ojf, O^i an <t an( i i s a l so red-shift dependent. 
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1. Introduction 

The concept of inertia has been one of the most debated topics of classical physics, 
starting with Newton's ideas of absolute space and that of inertia as an intrinsic 
property of matter devoid of any external influence. The notion of absolute space 
was criticized by Leibniz and later Bishop Berkeley who claimed that it is meta- 
physical. They were followed by Ernst Mach^who in 1872 rejected the existence of 
absolute space in favour of relative motion with respect to a "fixed" frame provided 
by the matter distribution in the universe, and claimed that it is the acceleration 
relative to this frame that determines the inertial properties of matter. This is the 
essence of Mach's Principle, a term coined by Albert Einstein in 1918, which says 
that "the inertial force which acts on an accelerating object is due to its interac- 
tion with all matter present in the rest of the universe" . There are indeed several 
interpretations of Mach's principle and arguments on whether the general relativity 
theory is Machian or not!^ - ^ Besides Mach's principle there are other approaches 
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towards the explanation for the orgin of inertia which describe inertia as a local 
instead of a global phenomenon. Among these one can mention the suggestion that 
inertial forces result from the interaction of matter with electromagnetic fluctua- 
tions of the zero point field, ^ (which is assumed to be homogeneous and isotropic, 
a property likely also related to the large scale distribution of matter in the universe 
and thus to the same frame envisioned by Mach) or the attribution of inertia to the 
interaction of a particle with its own field, namely the self-force.^ 

The implementation of Mach's principle as a quantitative law was done by D.W. 
Sciama ^ in 1953, by applying the formalism of electrodynamics to gravitation. 
According to Sciama the gravitational held of a moving particle in the universe can 
be expressed in terms of a "gravoelectric" part E and a "gravomagnetic" part B. 
These can be written in terms of the retarded scalar and vector potentials 4> and A 

by 

E = -grad<MI) l£ 

(1) 

B = curl A, 

where 4> is taken to be the Newtonian potential and A = ((/>/c)v(i); v(t) being the 
velocity of the particle. 

According to this law, the gravitational force between two masses mi and mi 
with relative acceleration a can be expressed by 

Gm 1 m 2 , Gm 1 m 2 a 

b 2 1 2 ' W 

r err 

in analogy with the Lorentz force, where the acceleration dependent term is referred 
to as "inertial induction" and is the result of acceleration in helds that obey retar- 
dation. Since the universe can be taken to be isotropic on the large scale, the 1 jr 2 
term will cancel when integrating over all masses in the universe to find the result- 
ing inertial force on an accelerating test particle of mass m. This force is therefore 
given by^ 

EGmMa Gmap f f f 1 , 
—2— = — ^ / / / Z&> (3) 



all m 



r 



where the region Rjj is taken by Berry (see Ref. [7]) to be the observable universe 
having density p given by that part of the universe which is contained within the 
Hubble sphere with radius c/H, representing the distance where galaxies are reced- 
ing with the speed of light. Hence 

_ ^Gmap [ C / H f2nGp\ 

Substituting values for p and H in the coefficient of ma in equation (j4|), Berry 
obtained a value of ^ instead of unity. He explained that this discrepancy results 
from the uncertainty in p which does not include dark matter and the increase in 
masses of distant galaxies by the relativistic factor (1 — v 2 /c 2 )~ 1/>2 . 
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An attempt to generalize Mach's principle to curved spaces, particularly the Fried- 
mann universes, was done by R. L. Signoref^HSl by using Sciama's law with a po- 
tential 



where R(t) is the scale factor, instead of the above Newtonian potential. Moreover 
instead of limiting the volume of the universe to the Hubble sphere, the above inte- 
gration is extended to the causal sphere of the accelerating test particle represented 
by the particle horizon of the particle. In the case of the spatially flat Einstein-de 
Sitter model this potential takes the simple value of — c 2 , and then ([1} gives the re- 
sult F = ma, in accordance with the equivalence principle. In the case of a spatially 
curved Friedman model the calculation leads to F = Xma where A > 1 is time de- 
pendent. This means that either the time dependent inertial mass is larger than the 
gravitational mass, or else, the gravitational constant G changes with time to keep 
the ratio of the gravitational to inertial forces constant. However Signore's work is 
based on two assumptions which in our opinion limit the validity of the obtained 
results. The first one involves the integration in ([5]), which as in Berry's calculation 
(j4]), assumes that the density of the universe p is constant over its entire history. 
Moreover Signore claimed that in a curved space the field acting on an accelerating 
test particle is still given by the non-covariant flat space expressions in ([I]). 

In this article we relax these two assumptions and use Sciama's approach to find 
the inertial force on an accelerating test particle arising from the inertial induction 
of all matter in the universe. As done by Signore we consider the particle horizon 
instead of the Hubble sphere to mark the region of the universe that influences our 
local dynamics!^ This represents the distance to the farthest objects that can be 
observed at any particular time, and the redshift becomes infinite for emitters on 
the particle horizon@ However unlike previous calculations by Berry and Signore we 
use a variable denisty p when integrating over the entire universe. Moreover instead 
of the expressions for the retarded potentials and fields in (JXJ) which are valid in 
flat space, we use general covariant formulations valid in a curved spacetime. The 
inertial force on an accelerating particle for the different types of Friedmann models 
is obtained in a form which is covariant with respect to coordinate transformations 
in the 3-space of these models. In section 2, we derive covariant expressions for the 
retarded potentials and fields for an accelerating test particle in the three cases of 
flat, open and closed Friedmann models. Then in section 3 we obtain expressions 
for the inertial forces and the results are summerized and discussed in section 4. 




(5) 



a We can receive photons emitted by receding emitters beyond the Hubble sphere («rec > c), 
because the Hubble sphere expands and overtakes these photons which end up inside the Hubble 
sphere where vrec < c. 
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2. Covariant Fields 



The study of electromagnetism in a general curved spacetime was done by DeWitt 
and BrehmeP] Later HobbJI^ obtained a formal solution for the electromagnetic 
potential^ = (0, A) in the case of a conformally fiat spacetime, and showed that 
in this case, the electromagnetic signals propagate only on the light cone without 
scattering off the Riemann tensor or radiation tail as in general curved spacetimes, 
thereby simplifying the solutions. 

Let us write the metric for the Friedmann universes in the form 

ds 2 = -dt 2 + R 2 {t)[dp 2 + fi(p)dn 2 }, (6) 

where 

{sinp, closed (k = 1) 
p, flat(fc = 0) (7) 
sinhp, open (k = —1) 

R(t) is the scale factor obtained by solving the Friedmann equations and dVl 2 = 
d9 2 + sin 2 8d(f> 2 . Then the potentials due to an accelerating point charge q moving 
along the path x' = z(t'), where (i',x') represent the retarded position, are given 

bjM 

R{t)a , , 

. q[G(*)d km , + H(*)n k n m ,]v m (i') {0) 

AkM = ' 



where 



F(tt) 



G(*) 



cot*, 


fc = 


1, 


1/*, 


k = 


0, 


coth*, 


k = 


-1, 


CSC *, 


k = 


1, 


1/*, 


k = 


0, 


csch*, 


k = 


-1, 


tan i*, 


k = 


1, 


o, 


k = 


o, 


tanh i*, 


k = 


-1, 



(9) 



and 



n k = R(t)V. x , ny =R(t')V; k ', (10) 
d km > = R{t)R{t')g hm ,, (11) 



b Greek indices take values from to 3; latin indices represent the spatial components 1 to 3. 
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cr = 1 + n k ,v k '; v k ' = dz k (t')/dt'. (12) 

The function \t(x, x') is the bi-scalar distanc^f) in 3-space measured along a geodesic 
joining the points x and x'. This can be expressed as L(x, x' , t) / R(t) , where L is 
the distance between x and x' at time t. The vectors n k and n k i in (fTU|) are unit 
tangent vectors at x and x' to the geodesic joining these two points, and the bi- 
vector parallel propagatoi^l d km i gives a relation between the components of a 
vector A 1 at x and the components of the same vector parallel transported along 
the geodesic joining x and x', 

A , — rl ,, A k 



A k — dkm' A 71 

The components of the electromagnetic field tensor 
b^ 



(13) 
A KV are given 



+ n k n ml )J r ' 



G 2 (*) / 



n k J 



")} 



(14) 



F km (x,t) = -^^y | G ^|^^ (dki>n m - d mV n k )J 1 ' 



where 



t' 



dtf' 



jl$pj{dki>n m - d mV n k )J 1 ' j ,(15) 
*(x,x'), (16) 



it R(t") 

expresses the relation between the retarded time t' and the time t; obtained by 
setting s — 0, and the four-current density is given by 

dz m ' B(t')v m ' 
J (x', t') = q% = q K > . (17) 
ds a 

The Lorentz force on a comoving test charge q having four- velocity u a in the Fried- 
mann models (|6]) , is given by the space part of the four vector 



i.e., 



/" = qF^u v , 



F = -qF l ° = qF°\ 



(18) 



(19) 



In order to obtain expressions for the fields in (1141) for different Friedmann universes 
represented by ([?])■ consider a point charge q at the retarded position (t', p') moving 



c In general two-point tensors, or bi-tensors, are used in the description of physical processes in 
which a cause at a point P 1 brings about an effect at a point P. Their indices (or arguments) refer 
to the points P 1 or P and are written with or without a prime over the index. For more details on 
bi-tensors see Synge (1960) in Ref. 1141 
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in the 6 = direction such that v' = (v cos9/(R(t'), -vsm6/R(t'), 0). To find 
the components of the parallel propagator dkm 1 in (HU) consider the change in the 
velocity vector as it is parallel transported radially along a geodesic in 3-spaces of 
constant curvature given in Then 

dv l = -T^dp. (20) 



2.1. Case k=0 

In this case the radial component of the velocity vector is unchanged under parallel 
transport and hence v 1 = v 1 so that g\y = 1. Also v 2 = —vsm8/R(t)p such that 
.922' = P and v 3 = v = 0. Taking the limit v « c the expression for the field in 
(fT4"]) with W = p becomes 

qR 2 (t') m , qp- k 

Fofe = ^T (3fem ' + P;fcP:m ' )a ~W (21) 

where a m are the components of the acceleration of the point charge at the retarded 
coordinates (t',p'). This expression reduces to — E given by (p} in the case of flat 
spacetime where R(t) — R(t') = 1. The total force on the accelerating point charge 
q from a charge distribution, can then be obtained by using (|19[) and integrating 
over all contributions. If the distribution is isotropic the 1/ p 2 term in (|21j) vanishes, 
and one only needs to consider the component of the force along the direction of 
motion, given by 



F = / p{x n )F m -^d 6 x n = - / p(x n )F 0t —d"x n , (22) 
J l a l J l a l 

where p(x n ) is the charge density of the distribution and 



2.2. Case k —1 

In this case v 1 — v 1 and v 2 = — v sin6 / R(t) sin p. Therefore gw = 1 and 322' = 
sin p. Also = p and the expression for the field (fT4|) in the limit v << c is given 

by 



qR 2 (t) , qp. k 

F 0k = -^T-- 5fem' + p;kp-,m')a - . 2 ■ 24 

R(t)smp R(t) sm p 



The component of the field in the direction of motion is then given by 



(25) 



19, 2013 9:21 WSPC/INSTRUCTION FILE text 



The Problem of Inertia in Friedmann Universes 7 



2.3. Case k =-1 



Again there is no change in the radial component of the velocity vector under 
parallel transport and so g\y = 1. The parallel transported angular component is 
given by v 2 = — v sin 0/R(t) sinhp, and so g22' = sinhp. The previous expressions 
for the held and its component along the direction of motion are given (for v < < c) 

by 

F ° k = pf ( !l (9km' + P;kP; m ')a m ' - B J P * 2 , (26) 
R{t)smhp R(t)smh p 



and 

F oijh= 5^m- asin ^- (27) 



|a| R 2 {t) sinhp r 



Sciama's law of inertial induction can be generalized to the curved spacetime of the 
Friedmann models by using (f2~Tj) . (fM|) and (|26p to represent the gravitational held 
of an accelerating particle and (|2"2")l for the gravitational force on this particle from 
an isotropic matter distribution instead of the previous expressions in |T]) and (J2|) 
which are only valid in hat space. This is done in the next section. 



3. Inertial Force 



Using the coordinate transformation r — Rofk(p), where Rq represents the scale 
factor at the present time t = to the Friedmann models in ([6]) can be written in the 
form 

dr 2 



ds 2 



-dt 2 +A 2 (t) 



1 — Kr 2 



r 2 dil 2 



(28) 



where A(t) is the dimensionless scale factor A(t) — R(t)/Ro and k = k/R^. In these 
coordinates the expressions for the component of the field along the direction of 
motion given by (|2"3"1) . (|2"5j) and (|2"7|> are given by 



a* _ qA{t') 
0i |a| A 2 {t)r 



a sin 



(29) 



for all models. Then using Sciama's law of inertial induction, the total inertial force 
on an accelerating particle of mass m in the Friedmann models is given by 

Gma f f f pA(t') 



F 



■ sin 9 y/g dr d9 d(f>, 



(30) 



where p is the density of the universe, and we set A(t) = 1 because the force on the 
accelerating particle due to matter at the retarded position (i',x') is calculated at 
the present time t = t$. Instead of the radial coordinate in ()30j) it is m ore convenient 
to use the redshift factor z = 1/A(t') — 1 which is related to r bjff^ 

-1 r 

=fk \ft&o~\ 



cm 



cO 



dz' 



E(z') 



(31) 
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where fk is given by ([7]) and Hq, Q c o represent the Hubble's constant and the 
current density parameter associated with spatial curvature respectively. Also for 
the metric (l28l) 



and 



_ A 3 (t')r 2 sm6 

V / g= n =2^; K 

VI — Kr 



Hq^cO 



1/2 



(32) 



E{z)= £fi i0 (l + z) n < . (33) 

_i(c) 

where rij depends on the source of the density parameter according to Table [1] 
Assuming that source of the energy density of the Friedmann models consists of 



Table 1. Values of rn for 
different sources in Eq. J33D 



Source 


Parameter 


m 


matter 


Qm 


3 


radiation 




4 


curvature 


n c 


2 


vacuum 


n A 






matter and radiation, the density p in ()30[) can be expressed in terms of the redshift 

by 




Substituting (f3T j) - (j34)) in (f30|) we get the inertial force on the accelerating particle 
for the three Friedmann models 

, n p r Mmo + Ml + z) [ f z dz> 1 

= 0; F = -ma / — / — — dz, 35) 

Jo {l + z)E(z) [Jo E(z')_ 

fc=1;F= -7ra/o (i+^w sm l^/o w)) d7 - (36) 

* = -l: F = -^L f ^o + »«,(! + *) sinh f ^ f (37 ) 

Unlike ((!]) where the integration is limited to the Hubble sphere, in our case the 
integration is extended over the causal sphere bounded by the particle horizon of 
the accelerating particle where z — oo. These includes the inertial contribution from 
all matter and radiation in the universe. 
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4. Results and Discussion 



In this note we have tried to compute the inertial mass of particles following the 
approach of Sciama, i.e. the analogy between electromagnetic and gravitational 
forces, assuming that the latter, just like the former are subject to retardation. 
As a result, we assumed that they yield, just like electromagnetism, a far-field 
force between material particles of masses mi, m<i proportional to Gmxm^aj c 2 r\ 
note that because of the 1/r (rather than the 1/r 2 appropriate for non-accelerating 
particles) force dependence, a spherically symmetric mass distribution yields a non- 
zero force at its center. Based on this analogy we followed Ref. [13] to compute the 
electric fields of accelerated charges in an arbitrary FRW space time and replaced 
the coupling constant by Gm\ to mimic the effects of gravitational force. Finally, by 
integrating this force over all matter in the universe we computed the gravitational 
force of the latter on a body of mass m in relative acceleration (with the universe!). 
This force is proportional to the acceleration a with a coefficient that depends on 
the cosmological parameters, in (philosophical) agreement with Mach's principle. 

Starting with the spatially flat Friedman model (k = 0), the inertial force on the 
accelerating particle given by (1351) in a matter dominated Einstein-de Sitter universe 
(f^Mo = 1: ^ao = Qro = 0) is given by F = — (l/3)ma. Although this coefficient of 
ma is larger than that obtained by Berry, it is still far from unity, and this implies 
that the total force on the particle cannot be explained only in terms of the inertial 
contribution from all the matter in the universe. In a radiation dominated universe 
(fiijo = 1, ^ao = &mo = 0) F = — (l/2)ma, while if we choose the currently favored 
cosmological parameters SImq = 0.3, f^Ao = 0.7 we get F — —0.23 ma. 

Having said this, we have to point out that the coefficient of ma in all three cases 
given by (|35[) - (|37l) depends on the current values of the cosmological parameters. 
Since these parameters depend on the redshift z according to the expression 

«i - r " jo(1 + z)Wi 1 (38) 



Ei( c ) n«(i + *)' 

where are given in Table [TJ this coefficient will be different at different epochs 
in the evolution of the universe. Therefore if Omo = 0.3, Qao — 0.7, then at z = 1 
^Mi = 0.77, Qai = 0.23 and the inertial force is F = — 0.31ma. This higher value 
of the inertial force is due to a higher percentage of matter at earlier times. For 
very large red-shifts the universe is radiation dominated with Qr — ¥ 1 and f^A 0, 
and the inertial force thus calculated becomes F — > (1/2) ma. It is not obvious at 
this point whether such a variation can lead to observable effects or inconsistencies 
with observations; we hope to be able to look at these issues in the future. 

For a spatially closed (k = 1) universe the effect of positive curvature leads, 
as expected, to a higher value for the inertial force on the accelerating particle. 
So for a slightly curved universe with Q c q = —0.1, f2j\/o = 0.3, J7ao = 0.8 we get 
F = — 0.24ma. As in the flat case the coefficient of ma is higher in the presence of 
radiation. So if Qmo = 0.3, S^ao = 0.7, O c o = —0.1, Qrq = 0.1, then F — — 0.33ma. 
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On the other hand in the spatially open (fc = —1) case, the negative spatial curvature 
reduces the inertial force. Taking Qmo — 0.2, Hao — 0.8, f2 c o = 0.1, we get F — 
— 0.19ma. The previous remark about the z-dependence of our results applies also 
to the curved Friedmann models. 

To summarize, in this article we have used the formalism of electrodynamics 
in curved spacetime to generalized Sciama's law of inertial induction to a curved 
spacetime represented by the Friedmann models, and obtained covariant expressions 
for the inertial force on a accelerating particle in these models. This extends previous 
work by Berry and later by Signore on the subject of inertia in Friedmann universes. 
As in Berry's case we have showed that the total inertial force F = —ma on an 
accelerating particle cannot be explained only in terms of inertial induction from 
the total matter and radiation in the current universe. We have also seen from 
([38| that the inertial force is redshift dependent and increases for earlier times 
until it reaches the asymptotic value of — (l/3)ma in the matter dominated case or 
— (l/2)mo in the radiation dominated case. 
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